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1. Introduction 

In quantum mechanics, the usual Rayleigh-Schrodinger Perturbation 
Theory (RSPT) is used for solving the eigenvalue problem of a Hamil- 
ton operator 

H = H + (3V (1) 

in terms of an unperturbed Hamiltonian Hq with known spectrum, a 
perturbation V and a coupling constant /3 > 0. The RSPT yields for 
the 7-th eigenvalue a formal power series 

(0) = + f3E[ I) + ... + 0»Ejp + ... (2) 

with real coefficients E„ , . Such series are often divergent, for instance 
for the quartic anharmonic oscillator (AHO) with Hamiltonian 

H = -^2+ x2 + ^ ( 3 ) 

on the real line. For this and closely related problems, there is a vast 
amount of literature since the quartic AHO may be considered as 
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the prototype of a zero-dimensional <3? 4 field theory, and thus, it has 
been used as test case for almost every important numerical method 
to solve the quantum mechanical eigenvalue problem. A recent ex- 
ample is the use of the generalized Bloch equation as an iterative 



method for the solution of the Schrddinger equation [18, 1£, pC| ]. For 
the quartic AHO, the n-th term a n in the RSPT series behaves as 
a n ~ - v / 247^(-3/3/2) ?1 r(n + 1/2) for large n §, £§ § and the series 
has zero radius of convergence in the variable (3. Thus, one has to sum 
such alternating divergent series. 

There are many summation methods that can be used in principle, 
a review of which is outside the scope of the present article. Consider 
a (formal) power series 



with partial sums 



f((3) = J2c j P j (4) 

j=0 



E C J P ■ (5) 



An important class of summation methods are nonlinear sequence trans- 
formations that transform the sequence {s n }^ =0 of partial sums to 
a new transformed sequence {s' n }^ =0 that is assumed to converge to 
the so-called antilimit of the divergent sequence {s n }. This antilimit is 
taken as the result of the summation process. 

We mention some examples of nonlinear sequence transformations. 
One is the famous epsilon algorithm |3l| It computes the upper half of 



the Pade table for the power series f(f3) according to p3| , 31] 

e$ = [n + k/k], {k>0,n>0). (6) 

Additionally, the epsilon algorithm is related to the Shanks transforma- 
tion |23| . Unfortunately, the epsilon algorithm is unable to sum several 
important perturbation series, for instance the RSPT series for the 



ground state of the octic anharmonic oscillator [ 29 1 . 

As a second nonlinear sequence transformation, we mention the S 
transformation of Weniger p8|, Sec. 8] that may be defined by the 



(7a) 
(7b) 
(7c) 
(7d) 



recursive scheme 
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5?((,{4K}) 



(7e) 



with ||, Sec. 8.3] 

q(k)_ (C + n + k)(( + n + k-l) 



((+n + 2k){( + n + 2k-l) 



(C>o) 



(8) 



The 5 transformation depends on an auxiliary sequence {w n }^Q of 
remainder estimates with uj n ^ 0. Several variants for the choice of the 
latter may be considered [l( 



13, 



For the choice ui r , 



»n+l 



originally proposed by Smith and Ford pq 
(or d variant in the notation of Weniger |2£ 

f 5W(c,K}) = 5( fc )(c,K},{ s 



one obtains the t variant 



i.e., 



71+1 



i}) 



(9) 



Eq. (|7|) is essentially one of the recursive schemes for the computation 
of the J transformation which is a rather general and well-studied 
sequence transformation that covers many of the most successful trans- 
formations as special cases |9, 1C, 11, 13, 15, 0]. The particular choice 
(H) corresponds to the fact that the <S transformation is identical to a 
special case of the J transformation, namely the case p = 3 of the P J 
transformation || |Tl|, |7j . 

Numerical comparison of results of the summation of the RSPT 
series for the ground state of the quartic AHO using the epsilon al- 
gorithm and the t-variant of the S transformation shows that the 5 
transformation is more powerful than the epsilon algorithm for the 
summation of this particular perturbation series that is rather strongly 
divergent. p9| 

An important additional consideration stems from regarding the per- 
turbation series as function for complex /3, i.e., one considers an analytic 
continuation of the formal power series to a function in the complex 
plane. Then, in principle, one has to allow for several branches of such 
a function. This question was studied intensively for the quartic AHO 
by Bender and Wu |3|, |2], |S| and Simon pi |5|. E(P) has a third-order 
branch point at (5 = that dominates the behavior for large /3, i.e., 
E(fl) = OiP 1 / 3 ) for -»■ oo. In the case of AHOs, often renormalized 
series are used that correspond — up to a factor — to a reexpansion 
of the original series in a transformed variable k. The relation between 
n and P for the quartic AHO may be obtained using Symanzik scaling 



x 



1 



a/4 



x plus a certain minimization criterion to be [^, 3C] 



P 



3(1 - k) 3 / 2 



(10) 
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After renormalization, one has to consider a Hamiltonian H(n) = 
Ho + kV with an unperturbed Hamiltonian Hq = —d 2 /dx 2 + x 2 and a 
perturbation V = (x 4 — 3x 2 )/3. Further, the eigenvalues of H(k) and 
H = H(/3) are related by E(k) = (1 — k) l l 2 E{f3). Besides making the 
renormalized expansion less divergent, the effect of such a renormal- 
ization is that the interval (0, oo) for (3 corresponds to a finite interval, 
e.g., to (0, 1) for k. 



Thus, E(/3) is a multiple-valued function for complex [24, 25|. This 
type of analytic structure is not easy to simulate using rational approx- 
imants in (3 that result on applying either the epsilon algorithm or the 
<S transformation (or other Levin- type sequence transformations Q) to 
the RSPT series in (3 of E(l,(3). Such a behavior of E(l,(3) can much 
easier be approximated by algebraic functions in f3 instead of rational 
functions. These algebraic functions are nothing but zeroes of polyno- 
mials with coefficients that depend on (3. For AHOs, this seemingly first 
was recognized by Cizek and coworkers who considered polynomials of 
a special structure, the so-called effective characteristic polynomials 
and applied these to sum the divergent AHO RSPT series Q . Effective 
characteristic polynomials have also been applied successfully for the 
extrapolation of quantum chemical many-body perturbation series [12] 



and been proven to be size-extensive, i.e., to scale correctly with the 
particle number |14| in this context. Later, Sergeev and Goodson [21] 



used more general polynomials that give rise to algebraic approximants 
|2], Sec. 8.6] for the summation of the ground state of AHOs. Algebraic 
approximants are generalizations of the quadratic approximants intro- 
duced by Shafer |22]| . These types of algebraic approximations will be 
discussed more thoroughly in later sections. 

All the methods mentioned above and most other methods for the 
summation of divergent perturbation series only use the partial sums 
or the coefficients of a single series, i.e., the RSPT series of one eigen- 
value for a particular value of the coupling constant (3. The rest of the 
methods may be classified as multipoint or multiseries methods. 

Multipoint methods are those where still a single eigenvalue E^> (/?) 
is approximated but several Taylor expansions at different points Pi 
are combined to construct the approximation. The most prominent 
example for such methods are multipoint Pade approximants H. It 
should be noted that the large-coupling regime (3 — > oo corresponds 
to k — ► 1 in the renormalized treatment and can be included in a 
multipoint treatment. For instance, Cizek and coworkers used two-point 
Pade approximants for this purpose ||. For details of such an approach, 
the reader is referred to the literature, e.g. Refs. | 
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In the present work, we consider multiseries methods. These are dif- 
ferent from multipoint methods since the aim is to approximate several 
functions depending on (3 simultaneously using their Taylor series at a 
common point, e.g., [3 = 0. This multiseries problem occurs naturally in 
quantum mechanical eigenvalue problems since the Hamilton operator 
normally possesses several discrete eigenvalues. Using the same splitting 
of the Hamiltonian in unperturbed operator and perturbation involving 
the same coupling constant j3 for all these eigenvalues considered as 
functions of j3, it is clear that the RSPT series of the form (||) for 
several values of / are not independent. In the context of rational 
approximation, this multiseries problem leads to simultaneous Pade 
approximants |2|, Sec. 8.1] or vector Pade approximants [§, Sec. 8.4]. 
On the other hand, combining the idea of a multiseries method and the 
approximation by algebraic functions, it is rather natural to construct 
a polynomial with coefficients that depend on j3 by using the coeffi- 
cients of several perturbation series simultaneously fl3|1 . The resulting 
algebraic approximations are called polynomial-type multiseries approx- 
imants and are the main topic of the present article. These are related 
to, but different from Hermite-Pade approximants as is explained later. 

The possible combination of the multipoint and the multiseries con- 
cepts leading to the construction of multipoint multiseries approxi- 
mants is not considered in the present work. 

The outline of the article is as follows. First, we will discuss the 
relation between effective characteristic polynomials and algebraic ap- 
proximants. Then, multiseries approximants are introduced and several 
special types of these approximants are defined. Some properties of 
multiseries approximants are derived. In the last section, numerical test 
results for the case of the unrenormalized quartic AHO are presented. 



2. Effective Characteristic Polynomials 

In this section, we sketch effective characteristic polynomials and some 
of their properties that are known in the literature, and point out 
that they are special cases of algebraic approximants. Also, the aim 
is to motivate the extension to a multiseries setting discussed in later 
sections. 

Consider a polynomial of degree N in E with coefficients that are 
polynomials in f3 of the form 

N N~j 

P N (E) =J2E j £ f NJ , k l3 k , f N , Ni0 = 1 . (11) 

j=0 k=0 
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Such a polynomial depends on (JV + 3) N/2 coefficients fN,j,k and is 
called an effective characteristic polynomial |j. 12, [14]]. The reason for 
the nomenclature is that characteristic polynomials of this form arise in 
the linear variational method in an orthonormal basis (f)j, j = 1, . . . , N 
for Hamiltonians H = Hq + (3V , and usually, the /'s are computed 
via matrix elements of H in this basis. For effective characteristic 
polynomials, however, the /'s are obtained from a perturbation series 

oo 

£(/?)= (12) 

by requiring that 

P N (E(fl) = 0(P N ^ N+3 ^ 2 ) (13) 

holds for /3 — » 0. This leads to a system of linear equations for the /'s 
with as many equations as there are unknowns. If this system possesses 
a solution, the /'s and thus, the effective characteristic polynomials are 
uniquely defined. For known /'s, the eigenvalues are approximated by 
zeroes of Pn that are denoted by IljVj, j = 1, ■ ■ ■ ,N and are called 
effective characteric polynomial approximants. The branch that repro- 
duces the input data is simply denoted by IIjv and called the physical 
branch. In order to discuss the dependence of the approximants on the 
coefficients Ej of the series ([H|), the explicit notation Unj [Eq, . . . , Em] 
is used where M = N(N + 3)/2 - 1. 

It should be noted that the physical branch can become complex. 
Then, for N = 2 both branches yield complex results, for instance. 
This, of course, is not reasonable for the computation of discrete eigen- 
values of the Hamiltonian that is a self-adjoint operator, and indicates 
a breakdown of the perturbative approach. Complex values of the 
approximants, however, make sense for the description of resonances 



or tunneling processes |l], Chap. 8], [17], and thus, there are problems 
where the possibility of obtaining complex approximations from real 
perturbation series is a desired feature. 

We now discuss some properties of these approximants. Note that 
equivalence and invariance properties of algebraic approximants are 
known from the literature [||, Sec. 8.6] and this is relevant since the 
latter approximants are generalizations of effective characteristic poly- 
nomial approximants as discussed below. 



As shown independently in the literature [12, 14] before the connec- 
tion to algebraic approximants was recognized, the effective character- 
istic polynomial approximant H2 is invariant under a repartioning of 
the Hamiltonian where Hq is replaced by the new unperturbed oper- 
ator (1 — o)Hq and a is some constant. Denoting the corresponding 



proceed_a.tex; 2/02/2008; 2:47; p. 6 



Extrapolation of Perturbation Series 



7 



coefficients of the RSPT eigenvalue series by Ej(a), one has 

U 2 [E ,...,E 4 } = U 2 [E (a),...,E 4 (a)} . (14) 

This means that the approximant is invariant under the repartioning as 
the true eigenvalue has to be. In Ref. [|l4|] it is shown that a particular 
effective characteristic polynomial approximant IIjvj for some N and 
j with 1 < j < N has the scaling property 

n jv, j [c Eq , . . . ,cE M ] = cIl N j[E , E M ] • (15) 

This is important to guarantee size-extensitivity, i.e., correct scaling 
with particle number in many-body perturbation theory fDf] . 

We now show that effective characteristic polynomial approximants 
are special algebraic approximants. These are generalizations of Pade 
approximants. Algebraic approximants are constructed via polynomials 
A® in (3 with deg(A^) = d k , k = 0, . . . , N such that for Q = -1 + 



N + J2k=o dk and for a given power series E(f3) as defined in Eq. ( |12| ) 
the relations 

N 

P N (E((3)) = Y,AW(/3)[E(f3)] k = 0(/3« +1 ) (16) 

fc=0 

are satisfied for (5 — > 0. Since this defines the polynomials only up 
to a common factor, we additionally demand A^ N \0) = 1 as a nor- 
malization. For given polynomials A^ k \ the algebraic approximants 
-^[^0,^1,...,^]^' are defined as the zeroes of the polynomial Pn according 
to 

Pn (E [do>dl ,..., dN]ij ) = , (j = 1, . . . , N) , (17) 

and depend on Q + 1 coefficients Eq, . . . , Eq that play the role of the 
input data. Again, one may define the physical branch as that one 
which reproduces the input data via a Taylor expansion in j3. 

Comparison with the definition of effective characteristic polyno- 
mial approximants reveals that the latter are nothing but the special 
algebraic approximants 

IIatj = -E'[Ar i Ar-i i Ar-2,...,0],j > U = 1) • • • ! N) ■ U 8 ) 

It should be noted that algebraic approximants 

— are special cases of Hermite-Pade approximants and are described 
by a quite elaborate mathematical theory Sec. 8.5, 8.6] 

— are useful for multi-valued functions (analytic continuation to an- 
other Riemann sheet) 
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have been applied successfully to anharmonic oscillators [21] as 
noted in the introduction. 



3. Polynomial-type Multiseries Approximants 

In this section, we define the polynomial-type multiseries approximants, 
and introduce three special cases (II-, D-, and P-type approximants). 
The relation to Hermite-Pade approximants is pointed out. An ex- 
ample is given for Il-type approximants. Further, it is proved that 
the polynomial coefficients of P-type approximants can be computed 
recursively. 

As indicated in the introduction, the key idea of the present work is 
to use coefficients of several perturbation series for different eigenvalues 
to compute the algebraic approximants. 

For the definition of polynomial-type multiseries approximants^, we 
consider a polynomial 

N 

P N (z) = Y,^ k HP)z k (19) 
fc=o 

where the coefficients are polynomials in /3 with deg(A^) = dk > 
for k = 1, . . . , N. They satisfy for (3 — > the order conditions 

P N (EW(J3)) = 0(/?^ +1 ) , (7 = 1, ... ,S) (20) 

for S given series E^ (/3) of the form (||) and S parameters Qi restricted 
by demanding 

S N 

£)(Qj + l) = # + (21) 

7=1 k=0 

The normalization condition is 4^(0) = 1. The latter and the order 
conditions (^) again lead to a system of linear equations for the coef- 
ficients of the polynomials with as many unknowns as equations 
if ( ^l|) is satisfied. Hence the polynomials A^ are uniquely defined if 
the linear system possesses a solution as is assumed in the following. 

Given the polynomials A^ k \ the multiseries approximants 
j4Qx,...,Qs] defined as the zeroes of P/v, i.e., 



Qs) 

■;d-N],j 



0. (22) 



These are called simply multiseries approximants in the following. 
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These approximants depend on (3 and the coefficients 

E o ,■■■■> Eq], I = 1,...,S. (23) 

We are especially interested in the case that S = N which means 
that as many perturbation series are used as there are roots of P/v, 
i.e., as there are approximants. In this case, the consistency condition 
becomes 

N N 

£Q/ = £4- (24) 

1=1 k=0 

Therefore, we define the following special cases: 
Il-type approximants with S = N 

Hi*] j = E $,Z-l$~l:$o]j > U = l, ■■■,*). (25) 



Here, Eq. (|24|) is obviously satisfied. 
D-type approximants with S = N for even N = 2D 

Here, Eq. (||) is satisfied: The left hand side is N(N + l)/2, the 
right hand side is D(N + 1), and both agree since N = 2D. 

P-type approximants with S = N 

P IQ] J = ^gfcSojj « U = 1,-..,N). (27) 



Here again, Eq. (yM is obviously satisfied. 



n-type approximants were introduced by the author [13], D- and P-type 
approximants are new. 

We note that these multiseries approximants 

— are related to, but different from Hermite-Pade approximants (2|, 
Sec. 8.5], To be more specific, we remark that the order conditions 
(|20| ) are different from those satisfied by the Hermite-Pade poly- 
nomials a^rnj (/5) that are contructed from the same S power series 
EW(P), 1=1,... ,S via, e.g., 

S Mi 

"iM + E E a Itmi {P)[E {I) m mi = 0{n (28) 

/=lm,=l 

for suitable Mj and r. This should be compared to Eqs. ([l9|) and 
(12T 
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— are useful for multi-valued functions (analytic continuation to an- 
other Riemann sheet), 

— have no convergence theory so far, and 

— are possibly cheaper than algebraic approximants since the nu- 
merical effort to compute several short perturbation expansions is 
much less than one long expansion. 

As a simple example, we consider the harmonic oscillator problem 
with Hamiltonian 

d 2 

H = --^ + (1 + (5)x 2 , H = H-(3x 2 (29) 
with eigenvalues 

EM = (27-l)vT + 0. (30) 

The coefficients of the perturbation series follow by Taylor expansion 
of the 

E W = 3 + ^-^ + l^-i|^ + ... (31b) 

Then, using the series for I = 1 up to E^ and / = 2 up to E^ as 
indicated by underlines in Eq. (31), we obtain the following results for 



n [2] ,i = 2 + p- 1^4 + 4/3 + 2/32, 

= l + i/3-i/? 2 + — /3 3 - — /3 4 + . .. (32a) 
2 1 8 M 16 128 1 K ' 

n [2]>2 = 2 + /3+i v /4 + 4/3 + 2/?2 

= 3+|^ + ^-^ + ^ + ... (32b) 

Here, we have underlined the terms that are correct in comparison with 
the exact results. Thus, we gain one order for E^ in comparison to 
the input data while there is no gain in E&. 

The following theorem shows that the coefficients in P-type polyno- 
mials leading to approximants -P[q+i] satisfy the equations for the poly- 
nomials leading to approximants P[q], and hence, these polynomials can 
be computed recursively. 
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THEOREM 1. Define 

N-l Q 

(e) = e n +j2 e1 J2 a k' Q) P k ( 33 ) 

for all Q and given N. The zeroes of the polynomial P^ are the P-type 
approximants P[qj , if for (3 — ► 

P ( n Q) (e ^ V) = 0(/? Q+1 ) , (I = 1,-..,N) (34) 

holds. Assume that for I = 1, . . . , N 

/Q+i \ 

P n +1) =O(/? 0+2 ) (35) 

holds for (3 — > whence the zeroes of this polynomial are the P-type 
approximants P[q +1 j. Define 

N-i Q 

(£) = + J2 EJ E 4'' Q+ V • (36) 

j=0 k=0 

Then for [3 —> 0, the relations 

p{Q) ^C^iVj =0(/5 Q+1 ) (37) 

AoW for 1=1,... ,N, and thus, P^ } = pj? ] . 
Proof. Put 

Q 

EV> Q \P) = ^#fy. (38) 

Then direct calculation shows 
p(Q) ( £(/,Q) (/3)) = _ J gO+i[ £ ;U.O)( i g)]JV-ijv£;gj. 1 

N-l Q 

j=0 k=0 
N-l 

3=0 

+0(f3 Q+2 ) 

= 0{f3 Q+1 ) (39) 
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Table I. Quadratic P-type Approximants for 
the quartic AHO with (3 = 1/100 



Q 


P [Ql,i 


P IQ],2 


2 


1 007375 


3 03646 


3 


1 0073736 


3 03653 


4 


1.00737368 


3.036525 


5 


1 007373671 


3 0365254 


6 


1.0073736722 


3.03652530 


7 


1.00737367206 


3.036525306 


8 


1.00737367208 


3.0365253043 


9 


1.007373672081 


3.0365253045 


10 


1.0073736720815 


3.03652530451 


11 


1.00737367208137 


3.036525304514 


12 


1.00737367208139 


3.0365253045131 


13 


1.00737367208138 


3.0365253045134 


14 


1.00737367208138 


3.03652530451334 


15 


1.00737367208138 


3.03652530451335 


oo 


1.00737367208138 


3.03652530451335 



□ 

Recursive algorithms are also known for certain algebraic approxi- 
mants |2lJ, and, more generally for the computation of Hermite-Pade 
polynomials ||, Sec. 8.5]. 



4. Numerical Tests 

In the test cases, we always treat the quartic AHO with Hamiltonian 
as defined in Eq. (||) for small (3 without any renormalization. The 
coefficients in all the perturbation series used below have been com- 
puted using RSPT up to order 20 in the wave functions, followed by 
Taylor expansion of the Rayleigh-Ritz expectation value with this wave 
function leading to the coefficients in the perturbation series for the 
energies up to order 41. The computations were done in MAPLE V™. 

In Table |[ we consider quadratic P-type approximants (N = S = 2) 
for various values of Q. The calculation of approximants P[q]j requires 
two Q-th order perturbation series for the two lowest eigenvalues. The 
results show that both eigenvalues are well approximated for Q = 15. 
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Table II. II- type Approximants for the quartic 
AHO with (5 = 1/100 



N 


n[jv],i 


H[AT],2 


2 


1.007371 


3.0376 


3 


1.0073738 


3.03650 


4 


1.007373667 


3.0365266 


5 


1.0073736724 


3.03652522 


6 


1.00737367206 


3.036525310 


7 


1.00737367208 


3.0365253040 


8 


1.0073736720812 


3.03652530456 


9 


1.00737367208140 


3.036525304509 


10 


1.00737367208138 


3.0365253045138 


oo 


1.00737367208138 


3.03652530451335 



Table III. D-type Approxi- 
mants for the quartic AHO 
with /3 = 2/10 



N 


D[N],1 


D[N],2 


2 


1.11 


3.8 


4 


1.117 


3.53 


6 


1.1181 


3.534 


8 


1.11826 


3.5377 


10 


1.11828 


3.5386 


12 


1.118291 


3.5389 


Exact 


1.118293 


3.5390 



The convergence of the approximants is somewhat faster for the lower 
eigenvalue. 

In Table ||, we consider II- type approximants for various values of N. 
The calculation of approximants n^jj requires N perturbation series 
with orders N, N — 1, . . . , 1. 

Only the results for the two lowest eigenvalues are displayed for com- 
parison reasons. The data show that II- type approximants converge 
somewhat faster than P-type approximants in the example treated. The 
results show that both eigenvalues are well approximated for iV = 10. 
The convergence of the approximants is somewhat faster for the lower 
eigenvalue. 
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For large j3, both II- and P-type approximants break down and start 
to produce complex approximants. 

In Table |TJ, we consider D-type approximants for various values of 
iV for a larger (3. The calculation of approximants -D[jv]j requires N 
perturbation series with orders N, N — 1, . . . , 1 (even N). Again, only 
the results for the two lowest eigenvalues are displayed for comparison 
reasons. The results indicate that this type of multiseries approximant 
can be useful for somewhat larger (3. Comparison values are taken from 
Ref. |1. 

In further studies it is planned to study these approximants also for 
renormalized perturbation series in order to see whether for these the 
range of applicability of the approximants is extended as is the case for 
several other methods. 



5. Summary 

— Effective characteristic polynomial approximants are special alge- 
braic approximants. 

— Polynomial- type multiseries approximants combine the informa- 
tion of several perturbation series and allow the summation of 
divergent series. They are different from Hermite-Pade approxi- 
mants. 

— Perturbation series for these multiseries approximants are less costly 
to calculate than for algebraic approximants (several short series 
vs. one long expansion). 

— The simultaneous calculation of several eigenvalues is possible. 
Higher eigenvalues converge slower. 

— The large number of variants of multiseries approximants should 
be explored further. 
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